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ABSTRACT 

We evaluate the branching ratios of the radiative decays B — > K*'y, and 
B — s> Ki'j in the framework of three-point function QCD sum rules. We find: 
T{B K*-f)/T{b S7) = 0.17 ± 0.05 and V{B Ki-f)/T{b 57) = 
0.30 ± 0.15. We also study, in the infinite quark mass limit, the connection 
between the form factor of the radiative B K*'y transition and the weak 
hadronic form factors entering the semileptonic B-meson decays. 
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Radiative S-meson decays of the type B K*j, with K* = K*{892), Ki{UOO), etc., 
are known to provide valuable information on the Standard Model (SM) at the quantum 
level In fact, these decays, induced by flavour changing h ^ s neutral currents, are 
controlled by the one-loop electromagnetic penguin operator which involves important 
SM parameters such as the top quark mass and the Cabibbo-Kobayashi-Maskawa matrix 
elements Vts and Vtb- At the quark level, it has been shown that QCD corrections 
lift the Glashow-Iliopoulos-Maiani suppression leading to a sizeable enhancement of the 
branching ratio BR{b 57). Indeed, preliminary experimental results have been recently 
reported from the CLEO II Collaboration 0: 

BR{B K*{892) 7) = (4.5 ± 1.5 ± 0.9) x 10"^ . (1) 



On the theoretical side, the calculation of the branching ratios of the exclusive pro- 
cesses induced by the b s transition is plagued by the usual uncertainties involved in 
determining weak hadronic form factors. Estimates are available e.g. from quark mod- 
els 1^, and from the Heavy Quark Effective Theory in the two extreme cases where the 
chiral symmetry is assumed for the light quark sector 0, or where the strange quark is 
considered to be heavy 0; the different estimates differ by up to an order of magnitude. 

Among the various theoretical approaches, QCD sum rules appear to be promising 
since they are, in principle, model independent. A calculation, in this framework, of 
B K*'y has been done some time ago |§] using the formalism of two-point functions. 
However, since B* cannot decay into B K* because of phase space, two-point function 
QCD sum rules provide information only on the coupling of the effective current to B*. 
Vector Meson Dominance must then be invoked, and an independent estimate of the 
strong coupling constant gB*BK* is required. The latter lies outside the framework of 
two-point function QCD sum rules and, therefore, it is model dependent. A possible way 
out of this model dependency is offered by the formalism of QCD sum rules for three-point 
functions. In this case the hadronic form factors entering the matrix element of S — *• K*j 
can be determined provided the leptonic decay constants of the B and K* mesons are 
known. The former can be estimated from two-point function QCD sum rules Q], as well 
as from lattice QCD |]10|, whereas the leptonic decay constant of K* can be extracted e.g. 



from the r-lepton decay data ITT . 



In this paper we pursue this three-point function QCD sum rule approach, which is 
similar in spirit to another calculation that can be found in the literature [0; however. 
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as it will be shown below, the numerical results of can be made more reliable, and a 
more accurate determination of the above mentioned decay channels can be obtained. In 
addition, we study here the connection which follows in the oo limit between the 

form factor oi B ^ K*'y and the semileptonic B-meson form factors. 

The electromagnetic penguin operator which governs the b —>■ s'j transition can be 
written, for << m^, as follows: 

Heff = Cm.e^'sa^^—^q^ , (2) 

where e and q are the photon polarization and momentum, respectively. The constant 
C contains the dependence on the Cabibbo-Kobayashi-Maskawa angles and on the heavy 
quark masses; neglecting rric it reads: 

C = ^4r.VZV,,F,{^) , (3) 



m 



w 



where the function F2 is given by [0]: 

^.W^,r'-{AW + i;^&.^]}. (4) 



In eq.(|) rj = as{fJ.) / dsimw) where /x ~ mb is the typical scale of the process, and F2{x) 
reads [g : 

F2[x} = - — ^ + -T-T7T-7T7 TtH^) ■ (5) 



x-l] 



3 12 12 2 x-l 



The function F2{x) depends weakly on the top quark mass, e.g. in the range 90 GeV < 
nit < 210 GeV, F2{x) increases from 0.55 to 0.68. Using = 4.6 GeV, tb = 1.4 ps 
and rrit = 120 GeV, this leads to the prediction of the inclusive b ^ s'j branching ratio: 
BR{b 57) = 2.2 ( I 1/0.042)2 ^ ^g^^ 

Let us now consider the decay B — > K*{892) 7. According to eq.(^) the amplitude for 
B{p) K*{p',r]) '~f{q,e) can be written as follows: 

A{B K*^) = {K*{p',v)\ Ji:^^ \B{p)) e, , (6) 

where the effective current J^^j is given by: 

J^ff = Cmbsa^''' ^^^^'\ ,b . (7) 



The two hadronic form factors Fi and G2 appearing in the matrix element in eq.(^, i.e. 

{K*{p',v)\saJ-^^^q''b\B{p)) = ze,.,^v*'fp"'F,{q') 

+ [r^;(ml-ml,)-ir^*.q)ip + p'),]G,iq') (8) 

{q = p — p') can be related using the identity cr'"'^ 75 = —^e^^^P" ap„. In fact, G2 = Fi/2, so 
that we need to compute only the form factor Fi at the kinematical point = 0. 

As usual in the QCD sum rules approach, we consider the three-point function corre- 
lator: 

T.,Mp\q) = (if I dx dy e'^^'-^-^yH0\TiUx)d,M4{y))\0) , (9) 

where the currents Ja and J5 have the same quantum numbers of the K* and B mesons, 
respectively, i.e. 

My) = q{y)il^h{y) , (10) 
whereas the operator 0^,^(0) is given by 

oA^) = m\(y,A^) , (11) 

( o'pu = |[7/x,7iy])- Using the decomposition 

Tc^^uQ" = ^e^^,^p'p'^T{p^, p'^, q^) (12) 

we can compute the invariant amplitude T in QCD through the Operator Product Expan- 
sion for p'^,p'^ large and spacelike Q obtaining a perturbative term and non-perturbative 
power corrections parameterized in terms of quark and gluon condensates. Keeping the 
terms up to D = 5, and in the limit = 0, we have: 



rpQCDi^^2 J2 „2 



(p , p , g 



0) 



3m^ 
8^ 



ds ds 



{s - s') 



p'^){s' — p'"^) 



< qq > 



[nit 



p'^){—p 



+ -^9 < I^Gq > 



mt 



+ 



mt 



2(m2 -p2)3(_p/2) 

1 



3(m^ — p2)2(— p'2)2 2(m^ — p2)2(— p'2) 



(13) 



^The extrapolation from large and spacelike to = is possible due to the large 6-quark mass. 



In eq.(|T^) we have not included the 0{as) correction in the perturbative term. As it 
will become clear shortly, the one loop perturbative contribution turns out to be rather 
small, on account of the restricted integration range in the dispersion relations. Hence, it 
is reasonable to expect the two-loop contribution to be negligible. Also, the D = 4 non- 
perturbative term from the gluon condensate is not yet known. However, we do not expect 
this term to influence much the results. For instance, in other QCD sum rule applications 
to heavy-light quark systems, where the Wilson coefficient of the gluon condensate has 
been calculated, it has been proven to contribute very little to the sum rules. Current 
uncertainties in the values of the leptonic decay constants will dominate the final error 
on Fi(0). 

The invariant amplitude T(p^, p'^, g^) can be related, by means of a double dispersion 
relation, to a hadronic spectral density which gets contributions from the lowest lying 
resonances, the B and K* mesons, plus higher resonances and a continuum: 

T^(/,P'',0) = —F^{0) 

nib (p"' — m^)(p'^ — m^.j 

+ higher resonances + continuum. (14) 

In eq.(0), fx* and Jb are the leptonic decay constants of the K* and B mesons, re- 
spectively, defined by: {0\ J^\K*{p',r])) = fK*mK*ri^, and {Q\Jr,\B{p)) = fBm^/mb. 
According to duality, the higher resonance states and the continuum contribution to the 
hadronic spectral function can be modelled by perturbative QCD. A prediction for -Fi(O) 
can then be obtained by equating the hadronic and the QCD sides of the sum rule in the 
duality region 

ml < s < So 
< s' < s = min(s — mf, s'q) , (15) 

where sq and s'q are effective thresholds. 

The convergence of the series involving the power corrections, and the dominance of 
the lowest lying states, can be improved by performing a double Borel transform, with 
the result: 
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f r r 

IB jK*mK* exp - 

irib 
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3m^ 
8^ 
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Fi(0) 
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— < gg > exp[- 



mt 



M2- 



1 — m 



2/ 



°Um4 ^ 3M2M'2 2M2 



(16) 



where we have adopted the notation g < qaGq >= itlq < qq >. 

Equation (|16D is the sum rule for -Fi(O). We use the following values of the parameters: 
< qq >= (—0.23 GeV)^, and itlq = 0.8 GeV^; the effective thresholds are varied in the 
range 33 — 36 GeV^ for s and 1.4 — 1.8 GeV^ for s'. The value of fx* can be obtained from 
the decay K*^^^: fx* = 0.22 ± 0.01 GeV; as for and we use mb = 4.6 GeV 

and, consequently, fs = 0.18 ± 0.01 GeV as computed in [jl4| 0. 



The standard numerical analysis of eq.(|T6]) consists in looking for a region in the 
(M2, M'2) space where the value of -Fi(O) is reasonably independent of the Borel variables 
M2,M'2, and of the thresholds sq, s'q. In addition, one may look for a hierarchy in the 
series of power corrections, although there are a few known examples where this hierarchy 
is absent, without affecting the reliability of the results. It turns out that, in the present 
case, the perturbative contribution in ([161) is smaller than the D = 3 non-perturbative 
term, due to the tiny integration region which corresponds to reasonable choices of the 
thresholds. For this reason we are forced to relax the criterion of the hierarchy among 
all the terms in the power series expansion. In this way, as shown in fig.(l), we find a 
stability region in correspondence to the value -Fi(O) = 0.35 ± 0.05 (the error is obtained 
by varying the parameters in their allowed intervals). This result allows us to estimate 
the fraction of inclusive 6 — >■ 57 decays represented by the exclusive channel B K*^: 

r = (— ) (1-^ Fi O p^ = 0.17 ±0.05 . 17 

i (0 ^ 57) ^rrib^ ^ ^ 



B 



Using the computed value of the inclusive branching ratio we find BR[B K*'~f) 
(4 ± 1) X 10^^, to be compared with eq.(|l|). 

*For a recent calculation of the beauty quark mass see ||l5| . 
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The result in eq.(|T^ is smaller than that obtained in |]T2[ by roughly a factor of 2. 



First, there is a technical difference between our determination of -Fi(O) and that of []T2|. 
While we start from the three-point function (P), which leads to a one-step determination 



of Fi{0) from the sum rule (0), the authors of |jT2| use the equations of motion to rewrite 
Fi(0) as a linear combination of two form factors which satisfy two different QCD sum 
rules. We find that this procedure is unnecessary, and that it may introduce artificial 



uncertainties in the results for -Fi(O). Second, the authors in [0 use — 130 MeV, a 



value which is not confirmed by recent lattice and QCD sum rules analyses both for finite 
and infinite heavy quark masses p, ^ . 



A calculation similar to the one described above can be done for the transition B 
Kij. Since the electromagnetic penguin is a spin- flip operator, the Ki meson is the ^Pi, 
l"*" orbital excitation in the kaon system and, in the framework of QCD sum rules, it can 
be interpolated by the current Ja = qial^s. 

As it is well known, there are two 1^ states in the strange particle spectrum, i.e. 
iri(1270) and fs:i(1400), whose full widths are: 90 ± 20 MeV and 174 ± 13 MeV, respec- 
tively The leptonic decay constants of these states can be estimated using experi- 
mental information on the hadronic decays of the r -lepton. From the measured branching 
ratio [|TT1 BR(r -> (is:i(1270) + 7^1(1400)) z/^) = 1.14 ± 0.5 x 10"^ and using |TT| for the 



fraction of the r decay into Ki(1270) and Ki{Um), the values: 0.48±0.38 and 0.66±0.35, 
respectively, we find: /a'i(127o) = 180 ± 81 MeV and /a'i(14oo) = 258 ± 89 MeV . These 
values are in agreement, within errors, with a previous estimate of the leptonic decay 
constant of the ^Pi octet jl^ , based on two-point function QCD sum rules. 



It is likely that the two observed states are a mixture of the ^Pi and ^Pi states 
induced e.g. by the K*Ti channel. Since the mixing angle is not well known, we assume 
that the ^Pi state is an intermediate state with mass 1300 MeV and leptonic decay 
constant = 200 MeV . Using the threshold Sq = 3 — 4 GeV"^ we obtain the value 
-Fi(O) = 0.5 ± 0.1, which corresponds to 

^^^,""^7^ =0.30 ±0.15. (18) 

This result indicates that this transition could be dominant with respect to the B K*'j 
decay. It is rather amusing to recall that while both the ^Pi and the ^Pi states are 
produced in r-lepton decays, as well as in hadronic reactions, the latter is forbidden in 
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the radiative decay of the B meson. Hence, the observation of the above decay channel 
could shed light on the mixing between the two states fri(1270) and i^'i(1400). 

Our analysis of the B K*'y transition can be concluded by discussing the limit 
rrib oo. Although this does not affect our estimates of the exclusive radiative decays, 
it is important in itself as it can provide information on the form factors entering the 
semileptonic decays of the B-meson. In fact, it is known that in the limit of an infinitely 
heavy b quark one can relate the hadronic matrix element of the effective operator re- 
sponsible for the transition 6 57 to the matrix elements of the weak currents between 
a heavy and a light meson |]18|, by using the flavour and spin symmetries of the Heavy 
Quark Effective Theory (HQET) W^. Indeed, one obtains the relation: 



= + n^') _ !!^£±^A,(g^)| (19) 

2'mB niB + rriK* 2mB J 

among the form factor Fi responsible for the transition B K*'~f and the semileptonic 
form factors V{q'^) and Ai{q'^) defined by: 



2V{q 



2\ 



{K*{p',r])\s^^{l-^,)b\B{p)) = ' e^uX.V*''p'p"' 

rriB + tuk* 

+ i{mB + mK*)Ai{q^)7]*^ + ... (20) 

{B = B^ or B^); in (|20| ) the ellipses denote terms proportional to (p + p% or q^. In using 
HQET to relate Fi to V and Ai one meets the following problem. The relation (|l^) is 
supposed to hold for ^ g^^^ = (m^ — mK*Y, i.e. at zero recoil where the predictions of 
HQET are in general reliable, whereas in i? K*'y we need Fi{q^) at q^ = 0. Moreover 
at ^ Qmax have the following scaling behaviour of the form factors: 

Fiiq'Lax) ~ V^b 

A.iqlJ ^ -L. (21) 
\/mb 



Therefore, the contribution from Ai in ([191) should, in principle, be neglected in the 
infinite quark mass limit, as it is next-to-leading. However, V^(g^) and Ai(g^) might have 
a different g^ behaviour and, hence, one cannot exclude that Ai(0) is comparable to V{0) 
in the m^, — > 00 limit. QCD sum rules can shed some light on this issue. 



The hmit nib — > c>o in the sum rule for Fi{0), eq.(|I^), can be taken after defining new 
low-energy variables y, yo and E as follows: 



s = ml + 2mhy 
So = ml + 2mhyo 

= 2mbE ; (22) 

and we recall that the scaling law for fs is given by: Jb = f/y/f^b, modulo logs. Before 
taking the infinite mass limit we observe that the term arising from the D = 5 condensates 
in eq.(|l^) seems to diverge. However, as discussed in pO| , ^ this divergence is only 



apparent. In fact the contributions from the D = 3 and D = 5 condensates should be 
resummed since they arise from non local operators: they are only the first two terms 
in a power series of the non local operator in the expansion parameter mg. However, 
the exact form of the non local operator is unknown, and the resummation procedure 
introduces an ambiguity. In |2y, Q an exponential form exp(— mgA) was adopted. 



but a behaviour such as (1 + ^^^)~" {n = 1,2,...) is also acceptable since the power 
series of both forms starts with the same term: 1 — m^A. We observe, though, that for 
B K*'y, and with mf, ^ 4.6 GeV ( no doubt a large mass), eq.(|l^) shows the feature 
that the non-perturbative contribution dominates over the perturbative one. Should we 
resum the non perturbative contributions in the exponential form, we would obtain the 
unlikely result that, when mf, oo, the dominant contribution to the sum rule would be 
provided by the perturbative term, contrary to the expectation based on the rather large 
value mb ~ 4.6 GeV. Therefore, we shall assume the form 

(1 + ^)-^ (23) 
n 

Moreover, if one adopts the reasonable criterion that when mb — > oo, the quark condensate 
contribution still dominates the sum rule, one gets n = 1 and therefore the sum rule takes 
the form: 



which shows that Fi[q ) behaves near the origin as follows: 



EM ^ -L . (25) 
^/rrib 
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The consequence is that the behaviour of Fi{q^) between ^ 0, eq.(^), and ~ Imax^ 
eq.(^), can be interpolated by the simple pole formula suggested by a dispersion relation 
dominated by the nearest 1~ resonance: 

FM') = , ■ (26) 

Let us now examine the validity of eq.([T|) at = 0. We can follow the same method 
discussed above to derive the behaviour of the sum rules for V{0) and Ai(0) in the m^, —>■ oo 
limit. We obtain: 

f^.m^,feM-j^mO) = -^^^l dyj^ ds s exp[-- - ^ 

171^1711, 

QEM^ 



b 

m 



m 



■K* 



M'2 



■)^i(0) 



ryo r+oo 

/ dy ds's' exp[— 
Jo Jo 



y_ 

E 



+ 



2 <^-'>(i-a9^r 



M'2- 



(2^ 



The behaviour in eqs.(p7|)-(|28|) is compatible with the scaling laws in eq.(|21|) only if 
the form factors V and Ai have different q^ dependence, for example is a simple pole 
analogous to ( P^D whereas Ai is nearly constant in q^. This is in agreement with the 
observation made in , where the q^ dependence of the form factors of the semileptonic 
B ^ p transition has been explicitly investigated. 



We notice, in closing, that using eqs. (p^ , |27| , pHf ), the relation ([TOD is verified at g = 



This is a direct confirmation of the conjecture put forward in [|T8| and in [gj], that in 
heavy-to-light transitions the heavy quark stays almost on its mass shell, and that the 
relations obtained at the zero recoil point remain applicable over the full q^ region. 
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Figure captions 



Fig. 1. Fi(0) versus the Borel parameters M^, M'" for sq = 36 GeV^, s'q = 1.8 GeV^. 
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